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Abstract 

We present a novel mechanism of phase-dependent electric transport in dif- 
fusive normal metal-superconductor structures. We provide a detailed theo- 
retical and numerical analysis of recent unexplained experiments essentially 
explaining them. 
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What is the resistance of a small normal structure adjacent to a superconductor? Su- 
perconductivity penetrates the structure provided it is short enough. A naive suggestion 
would be that the resistance vanishes. However, it is not so. The simplest way to see this 
is to relate the resistivity to the scattering in the structure. Normal electrons travers- 
ing the structure should undergo scattering even if their wave functions are distorted by 
superconductivity. 

If the structure is connected to two superconducting terminals having different phases, 
the resistance of the structure will depend on the phase difference. This provides the physical 
background for what is called Andreev interferometry. There is an outburst of interest in 
the topic. Different types of Andreev interferometers have been proposed theoretically [|2]|| 
and realized experimentally. [f§-|6| 

Andreev scattering reveals a significant difference between diffusive conductors and tun- 
nel junctions of the same resistance. Optimal interferometers are composed of tunnel junc- 
tions. H^j] Standard theory predicts that the zero-voltage, zero-temperature resistance of a 
diffusive conductor is not affected by penetrating superconductivity 0, although it is slightly 
modified at finite temperature, when the sample length becomes comparable with the su- 
perconducting correlation length in the normal metal, £ = Jd/ (7rT), D being diffusivity. 
Apart from this thermal effect, a small modification may arise from the weak localization 
correction. || 

However, the recent experiment || demonstrates a significant phase modulation of the 
resistivity in the diffusive regime at very low temperatures. The authors insist that their 
results can not be explained by existing theories. 

Below we present a novel mechanism that provides zero-temperature phase-dependent 
resistivity. Due to electron-electron interaction, a weak pair potential is induced in the 
normal metal, that leads to extra Andreev reflection occurring in the structure rather than 
in the superconductors. The relative change of the resistivity 5R/R ~ g, g being the 
interaction parameter, and may be of either sign depending on the sign of g. 

However, our numerical results shows that for the concrete structure its magnitude is too 
small to explain the data of Ref . . Careful analysis allowed us to conclude that actually 
the more trivial thermal effect has been observed. The data show excellent agreement with 
the results of our simulations. 

The most adequate theoretical description of the system is provided in the framework of 
the Keldysh Green's function technique elaborated in || for superconductivity. In the dif- 
fusive approximation, one first gets equations for the advanced (retarded) Green's function, 
which is a 2 x 2 matrix depending on coordinate and energy, G(x, e), with G 2 = 1, 

d x (DGd x )+i[H,G] = 0. (1) 

D being the diffusivity in the normal state, H = ea z + i(ReA(x)a y + lmA(x)a x ) .We will 
assume that the temperature is low enough and the size of the normal structure is large 
enough to satisfy T,D/L 2 -C As, A s being the energy gap in the superconductors. Then 
the boundary conditions for G take a simple form (a being Pauli matrices): G = o z in 
normal reservoirs, and G = o x sin + a y cos in a superconducting reservoir having phase 
0. 

The Green's function G determines the characteristics of the energy spectrum of the 
quasiparticles in the structure. To solve the transport problem, we must know how this 
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spectrum is filled by extra quasiparticles. The equation for the even-in-energy part of the 
filling factor reads 

d x {D(e, x)d x f(e, x)) + j(x)f(e) = 0. (2) 

The first term here describes diffusion of quasiparticles with a diffusion coefficient that is 
changed by the penetrating superconductivity, D(e) = DTr((G + G + ) 2 )/8. At zero tem- 
perature, G = G + and D remains unchanged. The second term describes absorption of 
quasiparticle excitations into the superconducting condensate, or, in other words, conver- 
sion of the normal current into the superconducting one. The coefficient 7 is proportional 
to the local value of the pair potential, j(x) = A(x)Tr[ia y (G(e, x) + G(— e, x))]. 

In a normal reservoir biased at voltage V with respect to superconductors /(e) = 
eV/AT cosh _2 (e/2T). This provides boundary conditions for (||). The current into a reservoir 
is determined by the local gradient of /. 

The common theoretical approach (see, for instance, |T0|,PJ3[]) disregards interactions in 
the normal metal, that leads to A, 7 = 0. Since for diffusive conductors at zero temperature 
the common non-interacting picture does not give the resistance change, we concentrate on 
the effect of A in the normal metal. This value can be calculated with 

A = g J detanh(e/2T)Tr[ia y (G A - G R )}/8. (3) 

This is the novel feature of the present approach. 

Let us first make a simple estimation of the magnitude of the effect. Since it is expected 
to be small, we solve Eq. ^| to first order in 7. This gives a relative resistance change 
5R/R ~ / ~fL 2 /D. At zero energy, 7 ~ A. In normal metal, the energies in the window 
~ D/L 2 contribute to A, therefore A ~ gD/L 2 . This results in a simple estimation for the 
resistance change, 

6R/R = gc($), (4) 

c being a dimensionless number depending on the geometry of the structure and on the 
distribution of the resistivity therein. It is important to note that c depends neither on the 
structure size nor on the absolute value of the resistivity, provided the temperature is low 
enough, T <C D/L 2 . The effect depends on the normal metal material by means of g and 
can be of either sign depending on the sign of g. If the geometry of the structure is well 
defined, the effect can be used for the direct measurement of interaction in normal metal. 

At qualitative level, the effect seems to explain the results of Ref . || . Indeed, the phase 
modulation of the resistance they observed at low temperatures was of the order of several 
percent and looked material-dependent, including the sign of the resistance change. This 
prompted us to detailed numerical calculations of the resistance of a concrete structure ( 
Fig. 1) which is very similar to the one used in Ref. ||. 

The structure consists of the current branch, the superconducting branch connected to 
superconducting reservoirs biased at the phases — $/2,$/2 respectively, and the extra branch 
made for technological reasons. The current flows as it is shown in Fig. 2 and the voltage 
difference between the points A and A' is measured. Owing to the symmetry of the struc- 
ture, the voltage and / distributions are antisymmetric with respect to the superconducting 
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branch, whereas A(x) and G(x) are symmetric. Superconductivity in the structure gets 
completely suppressed when <3> approaches n. 

One-dimensional differential equations (|l|), (§) shall be solved for each branch and then 
matched in the crossing points. Fist we calculate G(e) in all points of the structure. Due 
to the boundary conditions, it depends on $ in each point, thus providing the origin of the 
phase dependent effect. We obtain A(x) by integrating G over energy. Then we calculate 
7(0;) and make use of an analytical formula that relates the resistance change to j(x). Details 
of the calculations will be reported elsewhere. fl2| 

In Fig. 2 we have shown the calculated phase dependence normalized by its maximal 
value at 5(f) = 0, c(0)/c(O), contrasted with the experimental data for Ag. The phase 
dependencies look similar, but the magnitude of the effect cannot be satisfactory explained. 
According to the calculation, c(0) = 0.14. If we take an expected value g = 0.04 [|ll[] for 
silver, we would obtain 5R max /R = 0.003 whereas the experiment gives SR max /R = 0.1. If 
we do it the other way around and try to fit g from the experiment we end up with g ~ 0.7. 
That would bring silver to the rank of high temperature superconductors. 

This prompted us to check the possibility that had been rejected by the authors of Ref. . 
We have calculated the thermal effect on the resistivity neglecting interaction corrections. 
For a given geometry, the relative resistance change 5R/R is a function of L/C, and (p. Our 
results are presented in Fig. 3. As expected, the thermal effect vanishes both at low and 
high temperatures. The resistance at $ = reaches the minimum at L ss 3£. We have 
plotted the normalized phase dependence at L = 3£ in Fig. 4 along with experimental data 
and obtained a perfect match. The maximal values of the change are also very close to each 
other: (5R max / R)theor = —0.097 versus {5R max / R) exp = —0.11. From the estimations given 
in H we obtain indeed L/£ = 2.5 — 3 at T = 20mK. As we can see in Fig. 3, the thermal 
effect persists at rather high temperature, in agreement with the long high temperature tail 
observed in ||. 

All this allows us to conclude that the experiment can be perfectly described by 
existing semiclassical theory of superconducting proximity effect and thus to remedy the 
seeming discrepancy between theory and experiment. 

The remaining discrepancy for metallic samples can be easily understood if one take into 
account the sensitivity of the effect to a concrete geometry and to the inevitable inhomogene- 
ity of these ultra-small structures. This point of view is supported by large sample-to-sample 
fluctuations of the magnitude of the effect. |]] The understanding of the results for Sb sam- 
ples having high resistivity presents a certain difficulty. However, we notice that all essential 
features for these samples, such a small magnitude of the effect, sinusoidal phase dependence, 
positive sign, and the long high-temperature tail, can be well understood if the structures 
are not completely diffusive but contain tunnel junctions. || 

We are ready to present several conclusions. 

We have shown that the results of || can be perfectly explained within the existing theo- 
retical framework and be attributed to the thermal effect, provided a concrete experimental 
geometry is taken into account. 

We present a novel mechanism of phase-dependent resistance in hybrid normal metal- 
superconductor structures that works at zero temperature. The observation of this effect 
would allow a direct measurement of the interaction parameter in a normal metal. 
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Our results show that the observation of the weak localization correction || is a more 
difficult task than it was thought to be. At zero temperature, this correction will be masked 
by the novel effect we discussed, provided the interaction is not very small. At high tem- 
peratures, the correction would become comparable with the tail of the thermal effect at 
T > DRq/ L 2 R. For the structures used in Ref. || this would correspond to unreasonably 
high temperatures of 20K. 

The authors are indebted to V. Petrashov for numerous discussions of his results, B. Z. 
Spivak for his illuminating remarks concerning the interaction in normal metal, D. Esteve, 
M. Devoret, H. Pothier and S. Gueron for the discussion from which this work has emerged, 
A. V. Zaitsev and G. E. W. Bauer for helpful comments. This work is a part of the research 
program of the "Stichting voor Fundamenteel Onderzoek der Materie" (FOM), and we ac- 
knowledge the financial support from the "Nederlandse Organisatie voor Wetenschappelijk 
Onderzoek" (NWO). 
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FIGURES 



FIG. 1. The structure under consideration. 

FIG. 2. Normalized phase dependence of the novel effect. Squares correspond to the experi- 
mental data of Ref. 6. 

FIG. 3. Temperature and phase dependence of the thermal effect. The temperature is incorpo- 
rated into £ 2 = D/ttT. The phase difference changes from for the lowermost curve to 7r for the 
uppermost one with step tt/20. 

FIG. 4. Normalized phase dependence of the thermal effect at L = 3£. Squares correspond to 
the experimental data of Ref. 6. 
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